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Abstract 

In this study, the temperature dependencies of magnetic response functions 
of the anhydrous dihalides of iron-group elements are examined in the neighborhood 
of multi-critical points (tricritical, critical end point, double critical end point) and 
first order transition temperatures within molecular field approximation. Our find- 
ings reveal the fact that metamagnetic Ising system exhibits anomalies in the tem- 
perature dependence of the magnetic response functions for r < 0,3. In addition, 
we extensively investigated how an inter- and intra-layer exchange interaction ratio 
can influence magnetic response properties of these systems. Finally, a comparison 
is made with related works. 
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1 Introduction 



Intensive theoretical and experimental effort is devoted to investigate the 
multicritical phenomena for more than half a century. The tricritical point 
(TCP) is one of the first multicritical points investigated which can be roughly 
viewed as a point separating a second order transition line from a first-order 
transition line, at which the three coexisting phases simultaneously become 
critical. Itinerant ferromagnets [I], multicomponent fluid mixtures [2], penta- 
nary microemulsions [3], ammonium chloride [J], and 3 He — 4 He mixtures 
[5] are other systems that represent tricritical behavior. In addition, it is 
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shown that there exist tricritical points in an experimentally accessible three- 
dimensional space of the electric field, temperature and pressure in ferro- 
electrics [5J. On the other hand, a critical end point (CEP) appears when a 
line of second-order phase transitions terminates at a first-order phase bound- 
ary delimiting a new noncritical phase. At this multicritical point, a line of 
second-order phase transitions intersects and is truncated by a first-order 
phase boundary, beyond which a new noncritical phase is formed. Binary al- 
loys [7\, relaxor ferroelectrics [5], binary fluid mixtures [9J, ferromagnets [TU] . 
Ising random-field model [TT], and metamagnets 2()-"TI)1 are the physical sys- 
tems in which the CEP are common. In 1997, an extensive Monte Carlo (MC) 
simulation [13] presented the singular behavior on the first-order transition 
line close to CEP in a classical binary fluid [T^1I15|I16] . 

In addition to TCP and CEP, the double critical (bicritical) end point 
DCP appears where two critical lines end simultaneously at a first-order phase 
boundary. Double critical end points have been observed in binary and quasi- 
binary mixtures [T7j , and there is also some indication of a DCP in the meta- 
magnet FeBr 2 [T8fT9] . According to mean-field approximation (MFA), the 
next-nearest-neighbor Ising antiferromagnetic model, the layered metamagnet 
and the random-field Ising model present a DCP [2"0jllll21] . In addition, The 
MC simulations exhibited the decomposition of the TCP into a DCP and a 
CEP in three dimensional spin-1 Blume-Capel (BC) model [22] while in d = 2, 
only a fully stable TCP is observed |23j . Recently, Plascak and Landau studied 
a DCP in the two-dimensional spin-3/2 BC model via extensive MC simula- 
tions [21]. 

On the other hand, the behavior of the staggered and direct susceptibil- 
ities in the neighborhood of phase transitions has been a subject of experi- 
mental and theoretical research for quite a long time: In 1975, a two lattice 
model of antiferromagnetic phase transitions is discussed in detail, using the 
Gell-Mann-Low formulation of renormalization group methods and Wilson's 
e expansion. In this study, Alessandrini et. al. have obtained the disordering 
susceptibility and staggered susceptibility in terms of two-point function at 
zero magnetic fields and zero momentum [25]. Later, Landau has obtained 
Monte Carlo data for a simple cubic antiferromagnet with nearest- and next- 
nearest-neighbor interactions which reveal asymptotic tricritical behavior of 
the order parameter and high-temperature susceptibilities which are mean- 
field-like without corrections, in agreement with renormalization-group calcu- 
lations [26.27J. Using the high-temperature series expansion for the extended 
Hubbard model, Barkowiak et. al. have obtained the series to the sixth or- 
der for the staggered magnetic susceptibility and charge-ordered susceptibil- 
ity [2H]- Recently, Li et.al. studied the susceptibility of two-dimensional Ising 
model on a distorted Kagome lattice by means of exact solutions and the 
tensor renormalization-group (TRG) method [29J. In addition, magnetic be- 
haviors of — CU2V2O7 single crystals are investigated by means of magnetic 
susceptibility measurements [30]. Millis et.al. report measurements of the mag- 
netization and susceptibility of a series of samples of two different variants of 
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the molecular magnet Mnl2-ac: the usual, much studied form referred to as 
Mni2-ac and a new form abbreviated as Mni 2 -a,c-MeOH |31j . 

Metamagnetic materials are of great interest since it is possible to induce 
novel kinds of critical behavior by forcing competition between ferromagnetic 
and anti-ferromagnetic couplings existing in the system, in particular by ap- 
plying a external magnetic field. Magnetic materials that exhibit field-induced 
transitions can generally divided in two classes; (i) Highly anisotropic, (ii) 
Isotropic or weakly anisotropic. The phase transitions in anisotropic materials 
(class (i) ) are usually characterized by simple reveals of the spin directions 
which are in contrast with transitions in class (ii). The field-induced transi- 
tions in class (ii) materials are related to a rotation of the local spin directions 
[32] . Iron group dihalides; compounds such as FeCl 2 , FeBr 2 , FeCl 2 2H 2 0, 
FeMgBr 2 , CoCl 2 and NiCl 2 fall in the first class (SUES]- Some theoretical 
Hamiltonian models describing the behavior of iron group dihalides have been 
proposed. Monte-Carlo simulation [3637] and high-temperature series expan- 
sion calculations [38.39J have been performed on a simple cubic lattice Ising 
model with in-plane ferromagnetic coupling and antiferromagnetic coupling 
between adjacent planes (the meta model) and on the next-nearest-neighbor 
(nnn) model with antiferromagnetic nearest-neighbor (nn) and ferromagnetic 
next nearest-neighbor (nnn) interactions. Recently, a MC simulation has been 
performed on a quite realistic model of FeCl 2 in a magnetic field jH] and this 
typical metamagnet has also been treated by a high-density expansion method 
on a two-sublattice collinear Heisenberg Ising (s = 1) model with three- and 
four-ion anisotropy - 

Although much effort devoted the critical behavior of the metamagnetic 
systems, to the best of our knowledge, there has been no studies investigat- 
ing the temperature and field dependencies of direct magnetic and staggered 
magnetic susceptibilities in the neighborhood of multicritical critical points 
such as critical end point and double critical end point as well as first order 
transition points. 

The layout of this letter is as follows: The derivation of the expressions 
describing the mean field staggered magnetic and magnetic susceptibilities is 
represented in Section 1. The results describing the temperature and field de- 
pendencies of the direct and staggered magnetic response functions are given 
Section 3, and finally Section 4 contains the conclusions and discussions. 



2 Derivation of static staggered magnetic and magnetic suscepti- 
bilities of Spin-1/2 Metamagnetic Ising Model 

In order to obtain staggered magnetic susceptibility one should intro- 
duce a staggered external field H s to the system [20] whereas total (direct) 
magnetic susceptibility is the response of the total magnetization to a physical 
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external field H. Consequently the Hamiltonian of the spin-| metamagnetic 
Ising Model can be written as 



H = — JijSiSj — J SfeS 1 / — H^Si — H s ^2(—l) l Si 



(1) 



intra 



inter 



where the first sum refers to ferromagnetic couplings = J > between spins 
and Sj in the same x-y layers and the second sum denotes anti-ferromagnetic 
couplings J kl = J < between spins in adjacent layers. In addition, H and H s 
denotes the external physical magnetic field and external staggered magnetic 
field respectively. By making use of mean field approximation free energy per 
spin can be obtained as 

/ = |T[(1 + m a ) ln(l + m a ) + (1 - m„) ln(l - m a ) 



+ (1 + m b ) ln(l + m b ) + (1 - m b ) ln(l - m b ) - 41n(2)] 



(2) 



+\z 1 j'm a m h - -z 2 J{m\ + mf) - \H(m a + m b ) - \H S 



m h ). 



In the constant magnetic field distribution, the sublattice magnetization 
m a and m b are functions of the independent variables T, if, and H s so that free 
energy per spin represented by Eq.(|2]) is a non-equilibrium thermodynamic 
potential which depends on several order variables [2] . The equilibrium state 
corresponds to the minimum of / with respect to m a and m b . On the other 
hand, the conditions for a stationary value of / can be expressed as follows: 



= K{m a ) + ^(m a + m b ) - sZ# - ^H s = 



= K(m b ) + ^f(m a + m b ) - + ^H t 



and 



T 2 



d 2 f d 2 f 



( d 2 f y 

\dm a dm b J 



= k (m a )« (m b ) + t 1 k (m a ) + k (m b ) > 0, 
where ft(mj) and k (rrii) are 



Mm,- 



1 + 77lj 



^1^ /. . ^2^ \ 

^ (1 + ^7 K 



-«(-77li), 



(3) 



(4) 



(5) 



«'(m») 



T 



zi J' ' 



(6) 
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In order to investigate the behavior of the metamagnetic system in the 
neighborhood of phase transition points, it is more convenient to formulate 
the system in terms of total and staggered magnetization which are given as 
follows: 

m a + m b m a - m b 
m t = g ' ms = 2 ' 

Inserting Eq.Q in Eq.(|3]) one obtains the following mean field equations 
of state for the spin-1/2 metamagnetic Ising model on a cubic lattice as below 

sinh2(^) 



cosh 2(^==p) + cosh 2( 



knT 



sinh2( ga fc y ) 
cosh2(^g^) + cosh2(^gf^; 



Here, 5 = J — J' , t = J + J' , Z\ = 2 and z 2 = 4 . The metamagnetic spin- 
1/2 Ising model exhibits field-induced phase transitions. In Fig. 1(a) and (b) 
the variation of the antiferromagnetic and ferromagnetic order parameters are 
given in the field-temperature plane for a second order phase transition for 
r = 1.0. 

It has been shown in the extensive theoretical review by Kincaid and Co- 
hen that metamagnetic Ising model exhibits different types of phase bound- 
aries (20]. In this study, a Landau expansion of the free energy is performed 
and by a careful analysis of the signs of the coefficients, the possibility of dif- 
ferent phase diagrams has been revealed. Further, Moreira et. al. has extended 
this analysis considering terms up to twelfth order [43J. The Landau expansion 
consists in developing the mean-field free energy given by EqQ in a power 
series of the order parameter (m s ) which vanishes near the critical point: 

n 

y(T,H,m s ) = J2^k(T,H)mf (9) 

k=0 

According to the values and signs of the expansion coefficients one can distin- 
guish different types of phase transitions: (i) For t/> 2 = and ip4 < and ijj e > 
a first order transition appears, (ii) If ip2 — an d ^4 > then an ordinary 
critical point takes place, (iii) For ip 2 = and ^4 = and ^ > we expe- 
rience tricritical point as a function of the ratio of the exchange interactions 

(V = 5£) = 

trcp = (1 - 4)0, 
m TC p = y/l - trap, 
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where = z 2 J + Z\J and t t connects first and second order transition lines. 
On the other hand, for r = 0.3 and t = £tcp 4>2 — ?p4 — 4>6 = which 
denotes a higher order critical point (TCP) [20E3] . Consequently, topology 
of the metamagnetic Ising model phase diagram depends on the value of the 
ratio of the exchange interactions (77 = jry)'- (i) If r > 0.3, the transitions 
between the anti-ferromagnetic and paramagnetic phases are of first order 
at low temperatures and strong fields while it is of second order at higher 
temperatures. The two types of transitions are connected by a tricritical point, 
(ii) For < r < 0.3, the tricritical point decomposes into a critical end point 
(CEP) and a bicricital end point (BCP) with a line of first order transitions 
in between, separating two anti-ferromagnetic phases [2"0f4"0] . see Fig. 7(a)- (b). 
The staggered magnetic susceptibility of an metamagnetic sytem is 

din 



If one uses this definition and the equations of state given in Eq.(13), after 
some algebra the staggered magnetic can written as below 

,. C20i2 — CL22C1 /1 oN 

ii Sr ->U a 2 l a 12 — a 22 a ll 

Here, H r = H/J' ,H Sr = H s /J' , an, ai 2 , a 2 i, a 2 2, ci, and c 2 are given 
in Appendix A. Whereas direct magnetic susceptibility (x t) is the response 
function of a system to a physical field and magnetic susceptibility can be 
expressed as, 

dm t 

Xt = H l ™o3H r - (13) 
Following the similar steps we have used in obtaining staggered magnetic 
susceptibility, one obtains % t as 



r ^22^1 — bi 2 d 2 . , , 

Xt— hm {- — — - — \. (14) 



3 Results 



Fig.2(a)-(b) represents the behavior of staggered and direct magnetic sus- 
ceptibilities of spin- 1/2 metamagnetic Ising model in the neighborhood of 
TCP for r = 1.0. One can see from the figure that staggered susceptibility 
(Xs) increases rapidly with increasing temperature and diverges at the tricriti- 
cal point. Whereas there is a discontinuity in the direct magnetic susceptibility 
(Xt) at the TCP. At this point we should note that Zukovic et.al. has rep- 
resented a study on dilute metamagnetic Ising Model within effective field 
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theory [511553 which takes account the spin correlations. Comparing Fig. 12 of 
Ref. [54J with Fig. 1(b) of the present Letter, one can see that our results are 
in accordance with the results of effective field theory. It's important to note 
that, this behavior is in accordance with the existence of the discontinuity in 
Xt on the critical curve, confirming that the transition there is second-order in 
the mean-field approximation jjl] as well as effective field theory [54.55J. On 
the other hand, Fig. 3 illustrates the temperature variation of the tricritical 
direct magnetic susceptibility for various values of the ratio of the exchange 
interactions (r). One can see from this figure that the amplitude of the fer- 
romagnetic susceptibility rises considerably high values for r > 1.78. One of 
the characteristic behavior of the metamagnetic Ising model for strong anti- 
ferromagnetic case (r < 0.3) is the existence of the re-entrance phenomena. 
One can see this fact in the phase diagram of metamagnetic Ising model for 
r > 0.3. For high values of the magnetic field the system is in a disordered 
state for T r — > whereas there is a transition from disorder to order at a fi- 
nite temperature. In addition under goes another second order transition from 
ordered phase to disordered phase in high temperature regime ( see Fig. 7(a) 
and (b)). 

Fig.4(a) illustrates the temperature dependence of anti-ferromagnetic sus- 
ceptibility for H = H rDCP = 1.994 which corresponds to the double critical end 
point(DCP) of the spin-1/2 Ising model for r = 0.2. Here T rjVi at which stag- 
gered susceptibility diverges denotes the Neel temperature at which system 
undergoes a second order transition from paramagnetic phase to antiferro- 
magnetic phase. Whereas the second divergence of the staggered susceptibil- 
ity takes place at T rDCP . In addition one can clearly see that there exist an 
non-critical maximum at the ordered phase . This maximum corresponds to a 
anomaly in the multicritical behavior of iron group dihalides. 

It is important to emphasize that Selke has reported that there is two 
lines of anomalies in the field-temperature phase diagram of the spin-1/2 Ising 
model with in mean field theory at which which the temperature derivative 
of the total magnetization exhibits, at fixed field value a maximum below 
the transition point, as exemplified in Fig. 3 of Ref. [10]. In this study the 
anomalies are related to the competing ordering tendencies of the external 
field and the interlayer couplings in a metamagnetic crystal. We should also 
note that there are experimental data which emphasizes the anomalies for 
quite some time [20J. In addition, there have been various experimental study 
on the field-induced Griffiths phase in Ising type metamagnets such as FeBr 2 , 
FeCl 2 and Fe 1 _ x Zn x F 2 [oTIfST] . Fig.4(b) exhibits the temperature variation of 
ferromagnetic susceptibility for H = H rDCP . In this case the signature of the 
second order transition from paramagnetic phase to antiferromagnetic phase 
is a discontinuity in the direct magnetic susceptibility which is in accordance 
with the literature [H]. In addition, there exist a special multicritical point 
which separates the two different anti- ferromagnetic phases (AFI and AFII). 



7 



This special continuous phase transition is of fourth order and direct mag- 
netic susceptibility represents a discontinuity at the double critical end point. 
In addition direct susceptibility also represent a discontinuity at T rjV2 which 
corresponds to a regular critical point from antiferromagnetic phase to param- 
agnetic phase. Fig. 6 (a) shows the temperature dependencies of staggered and 
direct susceptibilities of the spin-1/2 metamagnetic system for r = 0.2 and 
H = H rcEP = 1.99176. At this value of the reduced physical magnetic field 
the system under goes two phase transitions of different character. The first 
transition is the CEP which is also of fourth order [20] and takes place between 
disordered phase at lower temperatures and anti ferromagnetic phase at higher 
temperature regime. One can easily observe from Fig. 6 that the staggered sus- 
ceptibility diverges at CEP whereas the direct magnetic susceptibility shows a 
discontinuity. Similar to the anomaly at H = H rDCP , both x t and x a makes 
non-critical maximums in the antiferromagnetic phase. In Fig. 6 (a) and (b) 
we have given the temperature variances of the magnetic response functions 
of the system for different constant reduced physical field values. One can 
see from these figures that the broad maximum in the ordered phase declines 
with decreasing the amplitude of the physical external magnetic field. Finally, 
the line of anomalies in the staggered and direct susceptibilities is depicted in 
Fig. 7. Here (T — H) x denotes the field and temperature values which both the 
staggered and direct susceptibilities exhibit a broad maximum in the ordered 
phase as exemplified in Figs. 4 and 5. Unlike the anomalies discussed by Selke 
the broad maximum does not diverges as one approaches the double critical 
endpoint. Further, the anomalies the magnetic response functions of the meta- 
magnetic Ising system disappears for the case r > 0.3 where the critical end 
point and the double critical end points emerges to a tricritical point. In this 
case there is no re-entrance in the phase diagram. 



4 Conclusions and Discussions 



In this paper the magnetic response of iron group dihalides in a field with 
weak ferromagnetic intralayer interactions and highly coordinated antiferro- 
magnetic interlayer couplings are studied within mean field approximation. 
The expressions that describe the staggered (anti- ferromagnetic) and direct 
(ferromagnetic) susceptibilities are derived by making use of mean field theory. 
The findings of this study can be summarized as follows: the direct suscep- 
tibility exhibits discontinuity not only at the second order transition point 
but also at multicritical points such as TCP, CEP, and DCP. In addition, the 
both magnetic response functions of the metamagnetic Ising model exhibits 
non-critical maximums in the ordered phase at the region of the H r — T r in 
where the system shows re-entrance phenomena. 
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6 Appendix A 



The coefficients an, a 12, ^21, a 22, Ci, and c 2 in Eq.(14) are denned as follows: 



an 



1 = 1 - (l - tank ( -2( mi -m s )+4r(m t+ m s )+H r +H ar ^ 2 \ ^ + fr^-l ^ 



012 = 1 - (l - ten/i ( -2(m t -m a )+4r(m t +m a )+if r +H Jr ^ 2 A ( _ 2 + 



fl 21 = -1 + (l - tank ^(^+ms+iri^-m s )+H r -H 3r y j ^ + 4^-^ 



(15) 



022 = 1 — ( 1 — ian/i 



-2(mt+m, s )+4r(m t — m s )+_g r — g s 



-2 + 4r)T" 1 , 



r = [ 1 - itm/l ^2(m t - ms )+4r(m t+ms )+H r+ g 3 ^ j y-1^ 



C2 = ( — 1 + tanh 



-2(mt+m a )+4r(mt-m 3 )+H r -H Sr \ 2 \ rp-l 
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7 Appendix B 



The coefficients &n, b 12 , b 2 ±, b 22 , d\, and d 2 in Eq.(14) are defined as follows 



6u = 1 - (l - tank ( -2(m t -m B )+4r(m t +m J )+H r +g at . ^ \ ^ + 



6 12 = 1 - ( 1 - tank ( - 2(m t -m,)+4r(m t +m B )+H r +g ar A ^ (_ 2 + ^T" 1 , 



621 = -1 + (l - tCWl/l ( -2(mt+m,)+4r(m t -m,)+H r -g 8t . ^ ^2 + 4r)T T 7 1 , 

(16) 



633 = 1 - (1 - t<mfc ( - 2(m t +m J )+4r(m t -m J )+H r -g 8r j 2 j (_ 2 + 4r ) T -l ; 
rfi = ^ _ ^ ^ (-2(m t -m 8 ) + 4r(m t+ m a )+g r+ H Jr ) ^ ^ 

d 2 = (l - tank ( (-2(m t+ m 8 ) + 4r(m t -m 8 )+g r -H Jr ) ^\ y-1 
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Fig. 2. Behaviors of staggered (a) and magnetic (b) susceptibilities each as a function 
of reduced temperature for H r = jr = 1.354 and r = J/J = 1.0. For these values, 
the system has a tricritical point (TCP). And also, staggered susceptibility increases 
rapidly with increasing temperature and diverges as the temperature approaches to 
the TCP. 
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Fig. 3. The temperature variation of the tricritical direct magnetic susceptibility for 
various values of the ratio of the exchange interactions. Here the arrows illustrate 
the phase transition temperatures. 




Fig. 4. The temperature dependencies of staggered and total susceptibilities in the 
neighborhood of double critical point (DCP) and second order phase transition 
point which takes place for the value of the reduced magnetic field H rDCP = 1.994 
for r = 0.2. Here the arrows illustrate the phase transition temperatures. 
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Fig. 5. The temperature dependencies of staggered and total susceptibilities in the 
neighborhood of the critical end point (CEP) and second order phase transition 
point which takes place for the value of the reduced magnetic field H rcEP = 1.99176 
for r = 0.2. Here the arrows illustrate the phase transition temperatures. 




Fig. 6. The behavior of (a) the staggered susceptibility \s (b) the direct magnetic 
susceptibility xt as a function of the reduced temperature, where T r = is for 
several values of reduced field, H r = 
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Fig. 7. The calculated mean field phase diagram of the metamagnetic Ising model 
for r = -y = 0.2 in the temperature- field plane, (b) Detailed phase diagram in 
the neighborhood of the critical end point (CEP) and the double critical end point 
(DCP). The dashed lines denote the anomalies in the staggered and direct suscep- 
tibilities. 
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